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The current status of the Standard Model predictions for the muon anomalous magnetic moment is described. 
Various contributions expected in the Standard Model are discussed. After the reevaluation of the leading-order 
hadronic term based on the new e^e~data, the theoretical prediction is more than three standard deviations 
lower than the experimental value. 



1. Introduction 



For a particle with spin s and magnetic moment 



2m 



(1) 



where e, m and g are the charge, mass and gy- 
romagnetic factor of the particle. In the Dirac 
theory of a charged pointlike spin-1/2 particle, 
(7 = 2, and QED effects slightly increase the g 
value. Conventionally, a quantity a = {g — 2)/2 is 
referred to as the anomalous magnetic moment. 

The electron and muon anomalous magnetic 
moments have been measured with a very high 
relative accuracy of 0.24 ppb [ IJ and 0.54 ppm [ 
[2], respectively. The theoretical prediction for 
is only mildly affected by strong and weak inter- 
actions providing a test of QED and giving the 
most precise value of the fine-structure constant 
a. In contrast, allows to test all sectors of the 
Standard Model since all of them contribute sig- 
nificantly to the total. 

Although the electron anomalous magnetic mo- 
ments is known much more precisely, a^^ is much 
more sensitive to new physics effects: the gain is 
usually ~ (m^/me)^ « 4.3 • lO''. The r lepton 
magnetic anomaly has even better potential, but 
because of the small lifetime of the r, it has not 
yet been measured with the best limits coming 
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from DELPHI [ 3 : -0.052 < < 0.013 at 95% 
confidence level. The sensitivity of the DELPHI 
measurement is still one order of magnitude worse 
than the predicted value of Or [ E] . 

Any significant difference of a°^P from a*j^ in- 
dicates new physics beyond the Standard Model 
(SM). It is conventional to write as 
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,liad 



(2) 



where the terms correspond to the contribu- 
tions of Quantum Electrodynamics (QED), elec- 
troweak (EW) and strong (hadronic) interactions. 
While discussing these terms and their precision, 
it is worth comparing them to the experimental 
result [|2]: 



(11659208.0 ±6.3) • 10" 



-10 



(3) 



The QED part is dominated by the lowest-order 
term, represented by one graph, first-order in a [ 
[5]. The number of diagrams for the second- and 
third-order terms is more than 100, but they (up 
to a^) are known analytically [ [SI [71 [H]- Tak- 
ing into account a recent more accurate numerical 
calculation of the terms [ [9] and the leading 
log terms[fl0 l [TT 1 [T2 l [T3] one obtains 

^QED ^ (116584718.09±0.14±0.04)-10-", (4) 

where the errors are due to the uncertainties of 
the 0{a^) term and a, respectively, and the value 
of a-^ = 137.035999084(51) from the latest mea- 
surement of Qe has been used [ [TJ [T3] . It is worth 
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noting that the 4-loop term equals 38.1 • 10^^° 
and is thus six times larger than the experimental 
uncertainty. Therefore, it is clear that its calcula- 
tion as well as that of the 5-loop one is necessary. 

The electroweak term is known rather accu- 
rately [ESI [16]: 



a 



,EW 



(15.4 ±0.1 ±0.2) • 10" 
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(5) 



where the first uncertainty is due to hadronic 
loops while the second one is caused by the er- 
rors of Mh, Mt and 3- loop effects. 

The hadronic contribution can also be written 



as a sum: 



had had,LO 



,had30 



,had.LBL 
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The dominant contribution comes from the 
leading-order term, which using dispersion rela- 
tions can be written as [ jT7l [T8] 



had,LO 



a 



where 



R{s) 



V 37r y 



ds 
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and the kernel K{s) grows from 0.63 at s = 4mJ 
to 1 at s — > oo, 1/s^ emphasizing the role of low 
energies. Particularly important is the reaction 
e^e~ — !■ TT+TT^ with a large cross section below 1 
GeV. Numerically, a)^'^'^'^ « 700 x lO"!" [[19], 
so we should know it to at least 1% to match the 
experimental accuracy. 



[26l [27[ [28] . The new p meson sample has an or- 
der of magnitude larger statistics and a system- 
atic error of 0.8%. SND measured F.,^ from 390 
to 970 MeV with a systematic error of 1.3%[[^. 
KLOE studied F-j^ using the method of radiative 
return or ISR[[30l[3Tl[32] at 590 < < 970 MeV 
with a sample of 1.5 • 10^ events and systematic 
error of 1.3%[ 33J. The ji^^j values from CMD- 
2 and SND are in good agreement. The KLOE 
data are consistent with them near the p meson 
peak, but exhibit a somewhat different energy de- 
pendence: they are higher to the left and lower 
to the right of the p meson peak. However, the 
contributions to from all three experiments are 
consistent. After further analysis of the data with 
a higher statistics (more than three million events 
selected), KLOE reports the jF^rl values closer to 
those of CMD-2 and SND and achieves a 0.9% 
systematic error [ jS^ . First results on the 27r 
contribution coming from BaBar will provide ad- 
ditional information on this channel [I35|. 

Also important is the contribution of hadronic 
continuum. Some idea about the contributions 
to (jhadxo coming from different energy ranges 
and their uncertainties is given in Fig. [TJ Below 




2. Evaluation of the hadronic term 

Several estimates of a}^'^'^^ appeared recently [ 
[20l [2T1 [22] based on the progress in the low en- 
ergy e+e~ annihilation and including the data not 
yet available previously [ [TH [23| [24] . 

As already mentioned, one of the largest contri- 
butions to a'l^'i'LO comes from the 27r final state 
(about 73%). Therefore, a high-precision mea- 
surement of the corresponding cross section is 
one of the main goals of low energy experiments. 
In addition to the previously published p meson 
data[[2S], CMD-2 reported their final results on 
the pion form factor F.,^ from 370 to 1380 MeV[ 



had,LO 



Figure 1. Various contributions to aj^' 



1.4 GeV practically all final states have been mea- 
sured with consistent results by the CMD-2 and 
SND groups in Novosibirsk [[36]. Figure [2] shows 
cross sections of various final states measured at 
CMD-2. Impressive results on various final states 
with more than two hadrons above 1 GeV were 
achieved by BaBar [ 37, 35] [Ml [101 [H] using the 
ISR method. In addition to measuring for the 
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Figure 2. Cross sections of various final states 
from CMD-2. 



first time the cross sections of new channels, they 
show that some of the previous results should be 
reconsidered. For example, the BaBar results on 
the process e+e~ 7r+7r~7r° [^7\ agree well with 
those of SND [|42] below 1.4 GeV, but are consid- 
erably higher and just inconsistent with those of 
DM2 [ 133] at higher energies. The overall con- 
tribution of the tt+tt^tt" final state from 1 to 
2 GeV (excluding the 7r"''7r~7r*' one) has been 
(2.45 ± 0.26 ± 0.03) • 10-^° before the BaBar data 
appeared and it becomes (3.25±0.09±0.01) -lO^^" 
if we replace the DM2 piece with a much more 
precise one from BaBar. 

Using the new data below 1.8 GeV discussed 
above in addition to the whole data set of [ 
[19l[23] for old experiments, one can reevaluate the 
leading-order hadronic contribution to a^. The 
data-based calculations [[^[5TJ[55] slightly differ 
by the integration method and the cut-off energy 
above which the predictions of perturbative QCD 
are used, but otherwise are essentially very sim- 
ilar. In Table [2] we show the results for different 
energy ranges following [ |^ . 

The theoretical error consists of 1.9 • 10~^° due 
to uncertainties of radiative corrections in old 
measurements and 0.7- lO^^*' related to using per- 
tubative QCD above 1.8 GeV. It can be seen that 
due to a higher accuracy of e+e" data the uncer- 
tainty of aJj^'^'LO is now 4.4 • IQ-'^^ (0.63%) com- 



pared to 15.3 • 10-1° of Ref.[[Ill and 7.2 • 10'^° 
of Ref.[[Mj- 

We move now to the higher-order hadronic con- 
tributions. Their most recent estimate performed 
in[|44] gives a{i=''^'"° = (-9.8±0.1) • lO^^" and has 
a negligible error compared to that of the leading- 
order one. 

The most difficult situation is with the light-by- 
light hadronic contribution, which is estimated 
only theoretically. Even the correct sign of this 
term was established quite recently [ i45j. The 
older predictions based on the chiral model and 
vector dominance [ 1461 I47j were compatible and 
much lower than that using short-distance QCD 
constraints [ HH] (see also[|3n])- Their approxi- 
mate averaging in[ [50] gives (jhad.LBL _ ^^2Q ± 
35) • 10~ii. Even higher uncertainty is listed in 
Ref. [ ini] who added some terms not taken into 
account in Ref.[[3H] to obtain (110 ± 40) • 10"". 
Two most recent updates give (116 ±40) • lO^^^ [ 
[52] and (105±26) • 10"" [|53]. It is very tempting 
to find an approach to estimate the light-by-light 
hadronic contribution from the data, like, e.g., 
it was done in Ref. [ [SJ , where CLEG single-tag 
measurements [ fSFj of 77* — > 7r°, 77, rj' were used to 
estimate the contribution from the pseudoscalar 
resonances. 

Using for the light-by-light term the result 
of Ref. [ [50] and adding all hadronic contribu- 
tions, we obtain aj^^^ = (693.1 ± 5.6) • 10"^". 
This result agrees with other estimations, e.g.,[ 
[H US ISS [53 [m [22] and its accuracy as well as 
that of the other recent data-based evaluations 
benefits from the new e"'"e~data. 

All separate contributions to a^^ are collected 
in Table [U Adding the QED, electroweak and 
hadronic terms, we arrive at the theoretical pre- 
diction of (11659180.3±5.6)-10-i°. The improved 
precision of the leading-order hadronic contribu- 
tion allows to confirm previously observed excess 
of the experimental value of over the SM pre- 
diction with a higher than before significance of 
3.3 standard deviations. Two other most recent 
evaluations also claim a large excess of 3.4<t [[21] 
and 3.1tT [ 122]. Results of the comparison are 
also shown in Fig. [3] For the first time during 
last years the accuracy of the SM prediction is 
slightly better than the experimental one. 
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Table 1 
Updated aj;^^, 



Ts, GeV 






504.6 ±3.1 ±1.0 




38.0 ±1.0 ±0.3 




35.7 ±0.8 ±0.2 


0.6-1.8 


54.2 ±1.9 ±0.4 


1.8-5.0 


41.1 ±0.6 ±0.0 




7.4 ±0.4 ±0.0 


> 5.0 


9.9 ±0.2 ±0.0 


Total 


690.9 ± 3.9exp ± 2.0th 



Table 2 

Experiment vs. Theory 



Contribution 




Experiment 


11659208.0 ±6.3 


QED 


11658471.81 ±0.02 


Electroweak 


15.4 ±0.1 ±0.2 


Hadronic 


693.1 ±5.6 


Theory 


11659180.3 ±5.6 


Exp. -Theory 


27.7 ±8.4 (3.3cr) 



How real is a very high accuracy of the leading- 
order hadronic contribution obtained above? We 
believe that we understand well the radiative cor- 
rections due to initial-state radiation and vacuum 
polarization, but should not forget that they arc 
numerically rather large and may reach ^ 20%, 
so their critical reanalysis and tests of the ex- 
isting Monte Carlo generators are needed. The 
situation with the radiative corrections due to fi- 
nal state radiation is not so well established, so 
we have to rely on the model of scalar electrody- 
namics and confront it with the data. This may 
increase the uncertainty of a)^'^'^'^ . There is also 
a question of double counting of the hadronic final 
states in the leading- and higher-order hadronic 
terms [[58]. 

One of the serious experimental questions is 
that of the missing states. An obvious candi- 
date is final states with neutral particles only. 



which have been badly measured before. Re- 
cent experiments in Novosibirsk in which the 
TT^^, ri'j , tt'^ TT^ 'J , riTT^-f final states were studied 
in the energy range from threshold to 1.4 GeV 
by CMD-2 and SND (see Rcfs.[ 59,i60] and refer- 
ences therein) showed that the cross sections are 
dominated by the p,uj,(t> mesons and thus the cor- 
responding contributions are properly taken into 
account. From the upper limits on nonresonant 
cross sections obtained in these papers we can es- 
timate that a possible, not yet accounted for con- 
tribution is a]','^'^'^'^ < 0.7 • 10-1". However, one 
should remember that there are no measurements 
at all of the cross sections of such channels above 
1.4 GeV although they are expected to be small. 

We have already mentioned serious progress 
with ISR studies from BaBar. The discus- 
sion of their effect on the Qfhad.LO estimation 
can be subdivided into two parts: new results 
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EJ,1995{e*e -based) 

210±16 

DEHZ,2003 (e*e -based) j 
180.9 ±8.0 

DEHZ,2003 (T-based) t 
180.9 ±8.0 

DEHZ,2006 (eV -based) , | 
180.3 ±5.6 ; 
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Figure 3. Comparison of from theory and ex- 
periment. 



on already measured states and studies of var- 
ious new final states. In the first part there 
are processes which cross sections are consis- 
tent with the older measurements and more 
precise, e.g., 27r+27r~, 7r+7r~27r°, .... There 
are also final states for which the cross sec- 
tions strongly differ from the older, less ac- 
curate measurements, e.g., tt~^tt~tt^, 6tt, .... 
In the second part there are final states, 
which have never been measured before, e.g., 
K+K-TT^n^, K+K-TT+TT-TT°, Att'^t], K+K-t]. 
Obviously, one should calculate what contribu- 
tion to a'l'^'^'LO comes from them and add it to 
the previous estimate. While doing that one 
should be very careful since any final state ob- 
served may be only a subset of more general pro- 
cesses. For example, the K~tt^tt~tt'^ final 
state may come from the process (fir], so that our 
estimate of the contribution to a]','^'^'^^ should be 
correspondingly divided by the relevant branch- 
ing fractions, in this case B{(l) K^K~)B{ri 
TT+TT^vr") = 0.1118, effectively increasing our es- 
timate of this contribution by a factor of 8.94! 
Fortunately, we are interested in exclusive cross 
sections only below 2 GeV and the new processes 
above usually have a rather small cross section in 
this energy range. The first estimate shows that 
these new contributions may increase ahad,LO 



(1~3)T0~^'^, only slightly decreasing the discrep- 
ancy between the theoretical expectation and the 
experimental result. 

In view of the new measurements of the cross 
sections of the processes with K^K~ and pi- 
ous in the final state one should carefully re- 
consider the contribution from the KKnir fi- 
nal states, which was previously estimated us- 
ing isospin relations [ [53]. Anyway, it is clear 
that we have to process new information thor- 
oughly and understand the size and accuracy of 
the continuum contribution below 2 GeV (now 
(62.4 ± 2.0 ± 0.5) • 10-^°) compared to that from 
the TTTT (now (504.6 ± 3.1 ± 1.0) • 10"^°). 

There is still no explanation for the observed 
discrepancy between the predictions based on r 
lepton and e"'"e~data[ '24]. For this reason we 
are not using r data in this update. One ex- 
pected that more light on the problem would be 
shed by the high-statistics measurement of the 
two-pion spectral function by Belle which prelim- 
inary results indicated to better agreement with 
e^e~data than before [[6T]. However, it turns out 
that while in a relatively small range of masses 
from 0.8 to 1.2 GeV the tttt spectral function mea- 
sured at Belle is below the ALEPH one, see Fig. 21 
this effect is compensated by the spectral function 
behavior at low and high masses, so that the re- 
sultant contribution to the hadronic part of the 
muon anomaly is about the same as before. On 
the other hand, a recent comprehensive analysis 
of the e+e^data below 1 GeV and those on the 2tt 
decay of the r lepton performed in Ref. [[62] shows 
that two data sets can be reconciled if mixing be- 
tween the p, u, (f) mesons is taken into account 
in a consistent way. 

Another interesting insight into the problem 
of the muon anomaly has been demonstrated in 
Ref. [[63] . The authors discuss the possibility that 
the observed discrepancy between the experimen- 
tal measurement of and its theoretical predic- 
tion may be due to hypothetical errors in the de- 
termination of the leading-order hadronic contri- 
bution. In particular, they show that if one tries 
to solve the problem by increasing the low en- 
ergy e^e~cross section by an amount necessary to 
bridge the discrepancy, this affects the hadronic 
contribution to the running fine-structure con- 
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Figure 4. Comparison of 27r spectral functions in 
T decays. 



stant decreasing the electroweak upper bound on 
the Higgs mass. As a result, it leaves a very nar- 
row window for the Higgs mass. They also showed 
that this scenario would require a serious revision 
of the e+e~cross section, which seems unlikely at 
the current level of experimental accuracy. 

What is the future of this SM test? From the 
experimental side there are suggestions to im- 
prove the accuracy by a factor of 2.5 at E969 
(BNL) or even by an order of magnitude at 
JPARC. It is clear that it will be extremely dif- 
ficult to improve significantly the existing accu- 
racy of the leading-order hadronic contribution by 
measuring the cross section of e+e" annihilation 
to better than 0.3% as required by future de- 
terminations of mentioned above. One can 
optimistically expect substantial progress from 
new high-statistics ISR measurements at KLOE, 
BaBar and Belle together with the more precise 
determination of R below 4-5 GeV from CLEO-c [ 
[64] and BES-III [|65] . Experiments are planned at 
the new machine VEPP-2000 now commissioning, 
which is a VEPP-2M upgrade up to ^5=2 GeV 



with L„iax = 10'^^ cm ^s ^, with two detectors 
(CMD-3 and SND) [ E^. A similar machine 
(DA$NE-II) is discussed in Frascati [H?]- We can 
estimate that by 2010 the accuracy of q}^^'^'^ will 
be improved from 4.4-10"^'' by a factor of about 2 
(to - 2.2 • 10~i") and the total error of 4.1 • 10-^° 
wiU be limited by the LBL term (3.5 • 10"^"), still 
higher than the expected 2.5 • 10"^° in E969. 

Let us hope that progress of theory will allow 
a calculation of a^^'^ from first principles (QCD, 
Lattice). One can mention here a new approach 
in the QCD instanton model[[68] or calculations 
on the lattice, where there are encouraging esti- 
mates of a)i^'i'LO, e.g.,[ :69i (667 ± 20) • 10-^° or 
attempts to estimate a^""^'^^^ [[70], see also Ref. [ 
[TT] discussing successes and difficulties of this ap- 
proach. 

In conclusion, I'd like to emphasize once again 
that BNL success stimulated significant progress 
of e~^e~ experiments and related theory. Im- 
provement of e+e" data (BaBar, BES, CMD-2, 
KLOE and SND) led to substantial decrease of 
the ajj^^^'"'-* uncertainty. For the first time the 
accuracy of the theoretical prediction is better 
than that of the experimental measurement. Fu- 
ture experiments as well as development of the- 
ory should clarify whether the observed differ- 
ence between a™P and a^^ is real and what conse- 
quences for the Standard Model and for possible 
New Physics [ 72 it implies. 
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